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Abstract: Aberrations limit the resolution, signal intensity and achievable 
imaging depth in microscopy. Coherence-gated wavefront sensing (CGWS) 
allows the fast measurement of aberrations in scattering samples and 
therefore the implementation of adaptive corrections. However, CGWS has 
been demonstrated so far only in weakly scattering samples. We designed a 
new CGWS scheme based on a Linnik interferometer and a SLED light 
source, which is able to compensate dispersion automatically and can be 
implemented on any microscope. In the highly scattering rat brain tissue, 
where multiply scattered photons falling within the temporal gate of the 
CGWS can no longer be neglected, we have measured known defocus and 
spherical aberrations up to a depth of 400 |a.m. 
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OCIS codes: (110.1080) Active or adaptive optics; (010.7350) Wave-front sensing; (110.0113) 
Imaging through turbid media. 
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1. Introduction 

For biological imaging applications, light microscopy always compromises between image 
quality, penetration and non-invasiveness. However, due to refractive index mismatch (e.g., 
between water and tissue for water immersion objectives) and inhomogeneity within the 
specimen, wavefront distortions limit signal and contrast in deep imaging, especially for 
confocal, multi-photon microscopy or full-field Optical Coherence Tomography (ff-OCT). 
With an Adaptive Optics system, the wavefront of illumination can be shaped to compensate 
the distortion so that a diffraction-limited focus is restored. However, a traditional Shack- 
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Hartmann sensor (SHS) can't be used easily in strongly scattering biological samples since 
there is no selection of the ballistic photons originating from the reference point at a given 
depth in the sample amongst all scattered photons coming from the bulk. 

To avoid this pitfall, several sensorless schemes have been implemented: (i) stochastic 
methods based upon genetic [1] or hill-climbing algorithms [2], in which strongly stained 
specimens and numerous iterations are needed; (ii) model-based approaches based upon 
modal wavefront sensing [3], in which the reduced number of measurements minimizes 
photobleaching and damage to the sample and which have been demonstrated in, e.g., 
confocal microscopy [4], structured illumination microscopy [5] and two-photon fluorescence 
microscopy [6]; (iii) pupil segmentation methods [7,8]; (iv) indirect aberration measurement 
through phase retrieval from images [9]. 

In contrast, coherence-gated wavefront sensing (CGWS) [10-12], which combined low- 
coherence phase shifting interfere merry (PSI) [13] and a real SHS or a virtual SHS (vSHS) 
[10-12], can measure the distortions of the wavefront in scattering samples and therefore 
allow their correction by adaptive optics. Aberrations can be corrected at a high temporal rate 
since the wavefront distortions can be obtained with a single measurement. These pioneering 
studies implemented CGWS in two-photon fluorescence microscopy, where a femtosecond 
pulsed laser was used as the low coherence light source and an appropriate amount of BK7 
glass was incorporated in the reference arm for the compensation of dispersion between the 
two arms. CGWS has been validated in a solution of scattering beads [10,12] and up to a 
depth of 200 |xm in the weakly scattering zebrafish forebrain [11]. Most recently, coherence 
gating with a real SHS has been demonstrated to reject unwanted back reflections [14] of lens 
surfaces far from the focus. A limit of this last approach is the difficulty to combine the 
sample aberrated wavefront and the reference arm wavefront. It prevents in particular using 
commercial SHS, since the collimated reference beam has to be recombined with the sample 
light behind the lenslet array in order to have interference to take place [14]. 

Here, we implement a new CGWS scheme based on a Linnik interferometer with a 
Superluminescent Light-Emitting Diode (SLED) as low-temporal-coherence light source. 
Compared to the previously described CGWS setup [10-12], its main advantages are the 
automatic compensation of dispersion between the two arms of the interferometer due to their 
symmetry and its possible implementation on any microscope due to its simple design 
[15,16]. In fresh thick highly scattering rat brain slices, for a medium range numerical- 
aperture (NA) and low-magnification (20 x /0.5) objective, we successfully measure up to a 
depth of about 400 \im a known defocus aberration, obtained by displacing the coherence gate 
(CG) position with respect to the actual focus position (AFP). With a high-NA and high- 
magnification (63x/0.9) objective, we measure up to a depth of about 200 |im the defocus as 
well as the 3rd order spherical aberration introduced by the CG displacement. In contrast with 
the previous publications using CGWS [10-12], where backscattered light was mostly 
ballistic, multiple scattering is found to be not negligible in the rat brain. In addition to the 
ballistic photons, the CG of the SLED source also selects photons that have experienced 
multiple scattering, if their time of flight lies within the temporal gate. The experimental 
results show that the multiply scattered photons selected by the CG begin to influence the 
wavefront measurements at shallow depths, of the order of the mean free path. However, we 
show that wavefront distortions can be measured at much larger depths, despite the fact that it 
relies on photons originating from a diffuse Coherent Volume (dCV), which is much larger 
than the Coherent Volume (CV) defined by ballistic photons only. CGWS measurement 
eventually fails, presumably when the dCV starts to deviate too much from the CV. 

List of abbreviations: NFP, nominal focus position; AFP, actual focus position; BFP, back 
focal plane; CG, coherence gate; CV, coherence volume; dCV, diffuse coherence volume; 
CGWS, coherence gated wavefront sensing; FOV, field of view; NA, numerical aperture; 
eNA, effective numerical aperture; OCT, optical coherence tomography; ff-OCT, full-field 
optical coherence tomography; PSF, point spread function; PSI, phase shifting interferometry; 
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SLED, superluminescent light emitting diode; SHS, Shack-Hartmann sensor; vSHS, virtual 
Shack-Hartman sensor. 



2. Materials and methods 



2.1 CGWS setup 

The experimental setup (Fig. 1) used to implement CGWS is a custom-made Linnik low- 
coherence interferometer (a Michelson interferometer with identical objectives in each arm) 
[15]. Two different pairs of water immersion objectives were investigated: medium range-NA 
and low-magnification (20 x /0.5) objectives and high-NA and high-magnification (63 x /0.9) 
objectives. The low-coherence source (SLED, center wavelength 750 nm, bandwidth 23 nm) 
is collimated, made circularly polarized with the polarizer and the quarter wave plate, 
expanded to overfill the pupil of the objectives and split with a beam splitter cube (BS1). The 
output of the interferometer is send to: (a) CCD2 conjugated with the object plane by LI, in 
order to position before each experiment the focus of the sample arm objective at the surface 
of the sample (the second surface of the coverslip), where the origin of the nominal focus 
position (NFP) is defined; (b) CCD1 conjugated with the pupil planes by the afocal telescope 
LI, L2; (c) a real SHS (not shown in Fig. 1) also conjugated with the pupil planes. The SHS is 
used for calibration and alignment of each arm individually while the pupil of the other arm is 
blocked. CCD1 is used to record the interference pattern of the light fields on the pupils pi 
and p2. Using PS1, the wavefront topography and hence the aberrations on the pupil pi of the 
sample arm are measured. A grade A cover slip (refractive index 1.51) is placed on the rat 




cover slip 



TSl 



manual translation stage 



Fig. 1 . Schematic of the experiment. Illumination: SLED (Superluminescent Light Emitting 
diode), Pol (polarizer), QWP (quarter wave plate), BS1 (beam splitter). Sample arm: Obj. 1 
(microscope objective), TSl (motorized linear translation stage allowing axial movement of 
the sample), coverslip (to protect the sample), 2D lateral manual translation stage. Reference 
arm: coverslip (to compensate dispersion from the coverslip in the sample arm), Obj. 2 
(microscope objective identical to Obj. 1), M (reference mirror), PZT (piezo actuator for 
phase-shifting), TS2 (motorized translation stage to adjust the reference arm length). Detection: 
LI and L2 (lens doublets), CCD1 and CCD2 (CCD cameras), BS1 (beam splitter), FD 
(adjustable field diaphragm of the microscope). The piezo, camera, TSl and TS2 are controlled 
by a PC running a modified version of the Light-CT software (LLTech). 
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brain slice to protect it. A similar cover slip is used in the reference arm to balance the 
dispersion and is slightly tilted to remove unwanted back reflections. A field diaphragm is 
placed between the two lenses of the afocal telescope to control and tune the microscope field 
of view (FOV). The FOV is calibrated for different values of the FD diameter both by using 
geometrical optics and by measuring it with an optical fiber translated in the focal plane of the 
objective. 

The coherence length l c of our SLED is 21.5 (xm, determined from its bandwidth Ak and 
central wavelength X as 0.88X7AL It corresponds to a CG length l c /2n of 7.9 (xm in the rat 
brain (refractive index n = 1.36 at 750 nm [17]). The CGWS mostly senses the light scattered 
within the coherence volume (CV), whose axial length (centered at CG) depends on the 
coherence length, which is defined by the FWHM of the self-coherence function [10]. 

The sample can be shifted with TS1 along the optical axis, bringing different sample 
layers into focus. At the surface of the sample, the focus position coincides with CG position, 
by alignment (Fig. 2(a)). When imaging into a sample whose refractive index n s is greater 
than the index «,- of the immersion medium, refraction at the surface causes the AFP of the 
objective to be shifted deeper into the sample with respect to the NFP (Fig. 2(b)). Because of 
the index mismatch between the refractive indices of water (used for immersion in the sample 
arm) and of the sample, the CV simultaneously moves in the opposite direction (Fig. 2(b)) 
[17]. By moving TS2 (Fig. 1), we can change the CV position and therefore either 
compensate this effect, or create a known additional defocus as well as some spherical 
aberrations (Fig. 2(c)) (see Appendix A). 

In our coordinate system, the AFP is taken as the origin and the positive direction points 
deeper into the sample. We define the CG position as the distance between the AFP and the 
CV center on the optical axis. Finally, as paraxial calculation is not valid for high-NA 
objectives, we locate the AFP with ray tracing method (see appendix A): a quasi point source 
of light is shifted along the optical axis and the AFP is defined as the position where the 
defocus aberration is zero (some higher order aberrations still remain). 




Fig. 2. Principle of the method used to create known aberration, (a) Initial state at the surface 
of the sample (no aberration), (b) When focusing into the sample, the index mismatch moves 
the CV and the AFP in opposite directions away from the NFP and introduces tractable 
aberrations, (c) By changing the reference arm length, the CV position can be further displaced 
to add or subtract aberrations. 



2.2 Wavefront reconstruction 

In order to extract the phase of the electric field corresponding to the CGWS-selected 
backscattered light on CCD1 (Fig. 1), we use four-step PSI [13,18] by shifting the mirror M 
mounted on the piezoelectric translation stage PZT (Fig. 1) in the reference arm. The mirror 
M is positioned at the focus of the objective during the alignment process using the SHS. The 
amplitude of the PZT is optimized to produce the optimal phase stepping of nil. Theoretically 
(see appendix A), the PZT movements may introduce some defocus and high-order 
aberrations, but only up to a RMS wavefront error less than 1 nm for both 20 x /0.5 and 
63x/0.9 objectives. We therefore omit this effect in our computation. The motion of the PZT 
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is synchronized to the frame rate of CCD 1. For each measurement, the corresponding electric 
field (defined as the "CGWS signal", and whose intensity is shown in Fig. 3(b)) is extracted 
from a recorded quadruplet of interferograms (Fig. 3(a)) by the four-step PSI algorithm. As 
the maximal frame rate of CCD 1 is 60 frames per second (fps), the phase can be extracted at 
the maximum speed of 15 fps. 



Fig. 3. Wavefront reconstruction procedure, (a) Raw image recorded on CCD1 (imaging the 
objective pupils), in a rat brain slice, for AFP located 160 um deep below the coverslip and CG 
= -15 um from AFP (20x/0.5 objective), (b) Corresponding amplitude of the electric field 
obtained by PSI. Bottom right corner, schematic of the virtual lens array, (c) Intensity 
distributions in the virtual image plane of the sublens in white of subfigure (b) obtained by 
discrete Fourier transform for 4 CGWS images obtained at 4 neighboring sample positions, (d) 
The wavefront reconstructed from the slopes of the centroids of the vSHS at a given position 
and then averaged over M = 5 neighboring positions (here corresponding to a Zernike defocus 
coefficient of 0.31 um). 

We use a vSHS [10,12] to unwrap the phase extracted by PSI and reconstruct the actual 
wavefront [15]. Though the speckle arising in scattering samples prevents the use in CGWS 
of traditional phase unwrapping methods, which are strongly limited by path inconsistencies 
due to singularities [19] and discontinuities [20], the vSHS overcomes these problems and 
allows reliable wavefront unwrapping in the presence of speckle [10-12,14]. The 
reconstruction of the wavefront was based on the algorithm implemented by Denk's group 
[10,12], in terms of the first 28 Zernike modes (up to 6th radial order) [21], which we briefly 
summarize here. The wavefront is numerically propagated through a virtual sublens array. For 
this purpose, the entire aperture is divided into a number of square-shaped sublenses (Fig. 
3(b)). For the 20 x /0.5 objective, we choose a sublens size of 20 by 20 pixels (49 sublenses 
across the pupil). For the 63 x /0.9 objective, the sublens size was chosen as 13 by 13 pixels 
(47 sublenses across the pupil). For each sublens, we perform a discrete Fourier 
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transformation of the electric field over its sub-aperture. With scattering samples, the 
amplitude distribution of the electric field on a sublens is speckled. As a consequence, the 
discrete Fourier transform also shows a speckle structure in the focal plane [1 1,22] (Fig. 3(c)). 
To reduce this speckle noise, we measure a number M of independent ensembles of scatters 
(Fig. 3(c)) [10,1 1]. In most experiments, we use M = 5 positions shifted on line with intervals 
of 3 (xm. In some experiments (see figure captions), a larger M value is used and the shifted 
positions are placed on a grid with 3 |im spacing. At each position, (i) a center of mass 
algorithm [23] is used to locate the centroid of the intensity distribution, which is preferable to 
other techniques [22], (ii) the slopes of the local wavefront are calculated, (iii) the wavefront 
(described with Zernike coefficients) is reconstructed by least-square fit to the array of local 
slopes, (iv) the obtained individual Zernike coefficients are finally averaged over the M 
neighboring positions. This spatial averaging reduces the maximal theoretical rate at which 
the phase can be extracted (typically to 2 to 3 Hz if the sample or the beam is moved fast and 
for M = 5). In our case, fast measurements were not required and therefore not performed (the 
sample was translated by a manual stage). 

3. Results 

3.1 Measuring known aberrations with CGWS in rat brain slices 

In order to show that CGWS is able to measure aberrations in vivo, it is necessary to generate 
well-defined aberrations. The rat brain itself has some aberrations (that remain to be 
measured), but in order to assess the quality of our measurement, we have chosen to impose 
known aberrations and check that CGWS is able to measure them correctly. Previous work 
has shown that moving the position of the CG relatively to the AFP induces known 
aberrations [24]. This relative movement can either be produced by changing the length of the 
reference arm (Fig. 2(c)) or simply by index mismatch when focusing in a tissue with 
refractive index different than water (Fig. 2(b)). The two main aberrations introduced are 
defocus and 3rd order spherical aberration and their magnitude can be calculated exactly from 
the objective NA (Appendix A). Here, the following experiment is performed: (a) change the 
relative distance between the AFP and the CG for different depths in the rat brain; (b) 
calculate theoretically the expected aberrations, which are depth-independent (see Appendix 
A); (c) measure the aberrations experimentally with CGWS and compare them to the 
theoretical predictions. Since the sample can itself create some unknown amount of 
aberration, we do not represent the absolute magnitude of the aberrations but their slope as a 
function of the relative position between the AFP and the CG. 

We observe that CGWS accurately measures the defocus slope at shallow depths for both 
objective pairs, since the measured slopes match the theoretical predictions (Figs. 4(a), 4(b), 
4(d) and 4(e)). When going deeper in the tissue, however, the measured slope is always 
smaller than the theoretical prediction, and rapidly diminishes (see, e.g., Figs. 4(d) and 4(e)). 
The critical depth at which the drop occurs depends on the objective pair used. While for the 
20x/0.5 objective we are able to measure the defocus down to 400 |xm with little or no 
deviation, the accuracy degrades very much more quickly with the 63 x /0.9 objective: at 300 
|im, we measure already only half of the expected slope and the measurement tends to a flat 
curve with a zero slope. 

We also compare the 3rd order spherical aberration measured by CGWS to predictions 
(Figs. 4(c), 4(f) and 4(g)). For the 20x/0.5 objective, the expected spherical aberration slope is 
too small to be accurately measured by CGWS (it corresponds to a RMS wavefront error of 
X/162 when the CG is displaced by 10 (xm). However, for the 63 x /0.9 objective, the 
theoretical 3rd order spherical aberration is much larger and CGWS is able to measure it 
accurately. Not only is the slope for shallow AFP depths close to the theoretical prediction, 
but also the drop in the measured slope occurs approximately for the same depth as for the 
defocus (compare Figs. 4(e) and 4(g)). 
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Finally, let us note that aberrations could not be measured at very shallow depths, for the 
AFP between 0 and 40 (im. Even though our CV has an axial FWHM of only 7.9 |itn, its 
shoulders pick up the strong reflection from the glass surface, which dominates the 
backscattering from the tissue at the center of the CV. 
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Fig. 4. Measuring known aberrations at different depth for the 20x/0.5 and the 63x/0.9 
objectives, (a) and (b): raw defocus measurement at different depths as a function of CG 
position for respectively the 20x/0.5 and 63x/0.9 objectives and theoretical curves, (c): raw 3rd 
order spherical aberration measurement for the 63x/0.9 objective at different depths as a 
function of CG position. In (a), (b) and (c), curves were vertically shifted for visibility, (d) and 
(e): from subfigures (a) and (b), defocus slope as a function of depth for respectively the 
20 x /0.5 and 63x/0.9 objectives, (f) and (g): slope of the 3rd order spherical aberration and 
comparison with theory for respectively the 20x/0.5 objective and 63x/0.9 objective. 

3.2 GCWS selects multiply scattered photons in rat brain slices 

In order to understand what limits the maximum depth at which CGWS can be performed, we 
need a better understanding of the selection of photons by the CG as a function of the depth in 
the tissue. To get an insight into which photons actually fall within the CG and which 
information they carry about the wavefront of interest, we measure the speckle size and the 
magnitude of the CGWS signal at different depths for the two pairs of objectives (Fig. 5). 

The speckle size, which is calculated as the FWHM of the autocovariance function of the 
electric field on the pupil [25], is the most important parameter to assess the size of the 
coherence-gated region. Assuming ballistic propagation to and from the CV and a random 
distribution of scatters within the CV, we expect the speckle size to depend on the CV 
extension. If the size of the CV is of the order of the extent of the point spread function (PSF) 
of the objective and if the CV lies around the AFP, scatters within CV cannot be resolved by 
the objective and produce therefore a wavefront without speckle. If the CV moves away from 
the AFP, or if the CV extends laterally (as when aberrations are present or when multiple 
scattering starts to be non-negligible), then we expect to see in the CGWS signal a speckle 
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structure, whose typical size is inversely proportional to the lateral extension of the CV. At 
the smallest depths, the size of the speckle is maximal when the defocus is minimal, as 
expected from the fact that the focus is the smallest at this position (Fig. 10 below). 

Moreover, the speckle size also decreases monotonically with the depth in the tissue (Fig. 
5(a), FOV 1000 |im and Fig. 5(b), FOV 530 (xm). This can be interpreted as the fact that the 
effective coherent volume, from which selected photons seem to originate, is enlarged in 
depth by multiply scattering. We denote in the following this volume the "diffusive 
Coherence Volume" (dCV), which presumably is larger laterally than the CV and has an axial 
extension larger than the coherence length. Monte Carlo simulation have explored the spatial 
extent of this volume [26], but as exact scattering properties of the brain (in particular the 
anisotropy of the scattering) are still ill known, it is hard to predict the shape of this volume. 
The lateral extension of the dCV (djcv) (estimated from the typical size of the speckle in the 
objective pupil (d spe ckie) by d dC v -1.22 Xf 0 bj/d spe ckie [27,28], where f obj is the focal length of the 
objective) increases monotonically and seems to saturate at large depths (Fig. 5(c), FOV 1000 
Urn and Fig. 5(d), FOV 530 |im). 

As the lateral extent of the dCV is several hundreds of |im at large depths, we checked 
whether the microscope FOV could limit the dCV and alter the CGWS measurement. When 
the CGWS system works, the FD (Fig. 1) is totally opened and it corresponds to the maximal 
microscope FOV (1000 (xm and 530 (xm for 20 x /0.5 and for 63 x /0.9 objectives respectively). 
By closing the FD, we measured the speckle size for several values of the FOV, calculated the 
corresponding dCV size (Figs. 5(a), 5(b), 5(c) and 5(d)) and computed the defocus slope to 
determine the FD range for which aberrations can still be estimated (Figs. 5(e) and 5(f)). As 
the FD is placed in the detection arm, it does not influence the geometry of photon trajectories 
entering in the sample. However, a photon exiting the sample with a large angle (as if arising 
from a point outside the FOV) is blocked by the FD. Multiply scattered photons selected by 
the CGWS gate can therefore be rejected by the FD. Closing the FD reduces the maximum 
dCV extension reached at large depth (Figs. 5(c) and 5(d)): some scattered photons in the 
dCV are blocked by the FD. This effect is present at large depths where multiply scattered 
photons contribute to CGWS signal. But perturbation of dCV by the FD is progressive and it 
can already be seen at depths smaller than the mean free path, (about 200 |xm in the cortex of 
adult rat in the far red [29-31]). The filtering of multiply scattered photons entering the 
CGWS gate by the FD is not precluding the correct estimation of the defocus slope (Figs. 5(e) 
and 5(f)), except for the very small range of FD size (FOV smaller than 180 (xm (20 x /0.5 
objective) or than 60 (xm (63 x /0.9 objective). Above these FOV, the maximal depth at which 
defocus slope is correctly measured is the same for all FOV. These curves show that the setup 
FOV doesn't determine the maximal depth of the CGWS measurement. 

The effect of the depth on the CGWS signal is also an interesting parameter. In classical 
ff-OCT, the drop in signal is expected to follow the exponential attenuation of the ballistic 
light from Beer-Lambert's law on a length scale given by half (because the light scattered at 
the focus travels twice the depth before exiting the sample) the mean free path (-100 |xm in 
rat brain). In Figs. 5(g) and 5(h), we represent the magnitude of CGWS signal as its total 
integrated intensity. The signal increases first until a depth of about 100 to 150 (xm and then 
decays slowly with a characteristic length scale larger than half the mean free path (700 (xm 
for 20x/0.5 objective, 250 (xm for the 63x/0.9 objective). This dependence of the CGWS 
signal on depth could be accounted for by the selection within the coherence gate of a large 
fraction of snake-like multiply scattered photons. This confirms that the CGWS signal is 
influenced by multiply scattered photons in rat brain tissue. 

A possible explanation of the failure of CGWS at large depth could be that the speckle 
structure of the wavefront could not be resolved anymore at large depths by the camera pixels. 
To address this issue, we analyzed the wavefront at the maximal depths of measurement for 
each objective while binning the pixels of the camera (Fig. 6) before using the vSHS. We 
observe that for both objectives the defocus slope is still reliably measured, as long as there is 
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Fig. 5. Transition from single scattering to multiple scattering in CGWS measurement within 
rat brain slices (male Wistar, 45 days old.) and influence of the microscope FOV. (a) and (b) 
Speckle size as a function of depth when CG is centered on AFP. (c) and (d) dCV lateral 
extension estimated from the speckle size as a function of depth when CG is centered on AFP. 
(e) and (f) Defocus slope as a function of depth, (g) and (h) The magnitude of CGWS signal as 
a function of the depth, (a), (c), (e) and (g): 20x/0.5 objective; (b), (d), (f) and (h) 63x/0.9 
objective. In each panel, the measurements are shown for four different diameters of the FD 
corresponding to four different FOV. 
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at least about 1.2 to 1.3 pixels per speckle on average. As the speckle size for both objectives 
is larger than 2 pixels at the critical depths where the wavefront measurement fails, this shows 
that the CGWS measurement is not limited by the camera sampling in our setup. 




equivalent speckle size (pixels) 



equivalent speckle size (pixels) 



Fig. 6. Influence of speckle sampling by the camera pixels on the CGWS measurement 
(wavefront obtained by averaging M = 40 positions). The slope of defocus is plotted at depths 
of 400 um (a) and 200 um (b) for respectively the 20x/0.5 and 63x/0.8 objectives as a function 
of the average number of pixels per grain of speckle, which was varied by binning the camera 
pixels before propagation through the vSHS (from left to right on the horizontal axis, the 
binning of the pixels is increased and therefore the speckle size in unit of pixels is decreased). 

4. Discussion 

Under the condition of single scattering regime, which is characterized by a ballistic transport 
to and from the CV and by a single backscattering event in CV, the wavefront distortions can 
be accurately measured by CGWS, as verified previously by [10-12]. We show in this paper, 
that for more strongly scattering structures such as the rat brain, multiple scattering starts 
playing a major role in CGWS. The photons selected by the temporal gate are not confined in 
the CV, but arise from a much larger volume referred to as the dCV. 

The existence of an extended volume selected by temporal gating explains the CGWS 
performance as a function of depth. Due to multiple scattering, extra photons in the dCV 
participate in the wavefront measurement and the deeper in the sample, the larger the dCV. 
This effect (a) decreases the size of the speckle in the pupil (making the wavefront harder to 
measure) and most importantly (b) perturbs the wavefront itself, since it is averaged over 
photons coming from a much larger volume than CV. 

Our main result is to demonstrate that the CGWS measurement of the known aberrations 
imposed in our experiment remains valid at large depths. CGWS could be performed indeed 
at depths much deeper than the depth at which only ballistic photons are selected by the CG. 
This result demonstrates the possibility to implement CGWS as wavefront sensor in a close 
loop system to improve microscopes working in scattering samples, such as two-photon 
microscopy or full-field OCT. 

However, at depths larger than a critical depth, the CGWS fidelity drops. We have shown 
that this critical depth was neither due to a spatial filtering of the dCV due to a limited FOV of 
the microscope nor due to a poor sampling of the small speckle grain by the camera pixels. 
Therefore the maximal depth is more probably related to the fact that above the critical depth 
the wavefront averaged over all the photons in the dCV does not carry the same information 
as the one originating from the ballistic CV. A better understanding of this limit and of the 
fact that it depends on the objective used will need a model of scattering in the tissue [26]. 

It is interesting to compare CGWS with ff-OCT or conventional OCT. Our CGWS setup is 
mainly different in two aspects: the spatial coherence of the source plays a major role here 
and the coherence gating is performed directly in the pupil of high-NA objectives. In regular 
ff-OCT with thermal light, photons are indeed selected in the coherence window both based 
on their temporal and on their spatial position: the multiply scattered photons are mostly 
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eliminated. In the CGWS setup, we measure the interference pattern with the reference in the 
pupil of the objective: we therefore select photons based only on their path length (the 
temporal coherence). As a consequence the multiply scattered photons that fall within this 
temporal gate participate in CGWS, even if they have been notably deviated. This results in a 
temporal gating only of the backscattered photons. We therefore select much more photons 
that just the ballistic ones. 

Finally, in the regime determined above, where known aberrations could be reliably 
measured (200 (xm depth for the 63x/0.9 objective, 400 \im depth for the 20x/0.5 objective), 
PSFs are calculated from the wavefront measured by CGWS using the first 28 Zernike 
coefficients whose tilt and defocus are numerically removed (Fig. 7). It is calculated at the 
experimental CG positions where the defocus is minimal (but never exactly zero). We observe 
that the PSF degrades notably over the range of reliability of the CGWS. Some aberrations 
(mostly defocus and some amount of 3rd order spherical aberration) stem from the coherence 
gating itself, since the CG position was only close to the optimal defocus position. Other 
aberrations could have been introduced by the sample itself (as, e.g., 3rd order spherical 
aberration due to index mismatch [17]). However, these results have to be taken cautiously, in 
particular we cannot say anything about the isoplanatism of the aberrations as long as the 
wavefront we measured is not applied to correct a full-field image. However they give an 
estimate of the degradation of the image expected when going in depth in the sample. 
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Fig. 7. Point Spread Function as extracted from the wavefront measured by CGWS for the 
minimal defocus position using the first 28 coefficients. The PSF was scaled over the whole 
gray values (0-255). We observe that the PSF degrades notably over the range of reliability of 
the CGWS. 

5. Conclusion 

We have implemented a new CGWS scheme based on a Linnik interferometer with a SLED 
as low temporal but high spatial coherence light source. Compared to the previous 
implementation of CGWS [10-12], its main technical advantages are the automatic 
compensation of dispersion between the two arms. Here the short coherence length for 
coherence gating is obtained from a broadband continuous source, with two symmetric arms 
to ensure that there is no phase delay between them for all wavelengths. Thanks to this simple 
design, a consequent advantage is his easy implementation on any microscope. Moreover, it 
offers in the future the possibility to modify easily the CGWS setup in terms, e.g., of its 
central wavelength or its coherence length. 

In fresh rat brain slices, we successfully measured up to a depth of about 400 (xm for a 
20x/0.5 objective and 200 |im for a 63^/0.9 objective known aberrations, obtained by 
displacing the CG position with respect to the AFP in the sample. However, measurement of 
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the speckle size and the CGWS signal as a function of depth in the sample demonstrates that 
the CG was not successfully rejecting multiply scattered photons even at shallow depths. This 
was attributed to a large amount of multiply scattered photons, which could have similar time 
of flight in the sample as the photons of the reference arm. CGWS could be directly 
applicable at shallow depth or in thin slices, where regular wavefront sensing methods fail. It 
would allow the implementation of adaptive corrections. At larger depth, our results show that 
CGWS allows the quantitative measurement of known aberration despite the selection of a 
large amount of multiply scattered photons. Its benefits for close-loop adaptive optics have to 
be demonstrated by coupling it to any microscopy modality and by measuring the image 
improvement. Finally, imaging at large depth may require the improvement of the rejection of 
the multiply scattered photons, e.g., by the increase of the scattering length, e.g., by the use of 
higher-wavelength light [32,33]. 

Appendix A: Analytical computation under sine condition of the aberrations for a 
refractive-index-mismatch sample covered by a coverslip 

For an objective operating under the sine-condition, an axial shift of the diffraction-limited 
focus corresponds to a defocus and to all orders of the spherical aberrations [34]. Binding et 
al. have proposed a concise formula to predict the aberrations for the mismatched index 
situation taking into account the spread between the NFP and AFP [17]. Here, we derive the 
analytical aberrations when a cover slip is inserted between the sample and the immersion 
water, as it is the case in most in vivo recordings when a glass coverslip is used to stabilize the 
brain after exposing it through a craniotomy. 

Under the sine condition, the principal surface of the objective corresponds to a sphere 
segment of radius nj (Fig. 8) [35]. To determine the wavefront from a point A (i.e., the single 
pass aberrations) on the optical axis, we imagine a point source at this position and trace its 
rays into the back focal plane (BFP) of the objective. Ray tracing restricted to the meridional 
plane is sufficient because of rotational symmetry. 
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Fig. 8. Geometry of schematic for ray tracing. Refractive index of water, glass slip, and sample 
are n b n s , n, respectively. If the sample was a pure water solution and if there was no coverslip, 
the focus in this aberration free case would be located at the position N on the optical axis. The 
origin 0 is located at the cross point of the optical axis to the second surface of the glass 
coverslip, the Z axis is defined along the optical axis OA with the positive direction pointing 
towards the sample (away from objective). The distance ON is noted d, the thickness of the 
glass coverslip T, the actual point source is A and the distance AN is noted Az. 
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Emanating from A at an angle a s with respect to OA, a ray is refracted and crosses the 
coverslip surfaces at P and B. The refraction angles in the glass coverslip and in the 
immersion water are a g and a, respectively. It crosses the BFP at E (at a distance r of the 
optical axis OZ). An auxiliary line BC is drawn perpendicular to the ray NC emanating from 
the NFP (N) at an angle a, with respect to OA. According to Fermat's Principle, all rays of a 
planar wavefront transmitted through the objective have the same optical path length when 
they reach their common focusing point in the BFP. Thus all the rays at an angle a, with 
respect to OZ in water arrive at point E, which is at a distance r=ntf sina, of OZ on the BFP. 

For the light emanating from the NFP N in the aberration free case, a spherical wavefront 
will be produced at the principal surface, which corresponds to a planar one at the BFP. The 
optical path difference W(f) between the rays originating from A (in the presence of sample 
and coverslip) and from N (in the aberration-free case, i.e., in the absence of sample and 
coverslip, with the full light path in water) and reaching the BFP in E, is equal to the optical 
path from A to B minus the optical path from N to C. 

The optical path from A to B is 

OP(A) = n s AP + n g PB = n s (d + Az)/ cos a s +n g T I cos a g . (A. 1) 
The optical path from N to C is 

OP(N) = n t NC = n, [(d + T) I cos a, + BD sin a t ] , (A.2) 
where BD = (d + Az) tan a s + T tan a g -(d + T) tan or, . Thus W(r) is 

W(r) = (d + Az)(n 2 -N 2 y 12 + T(n 2 g -N 2 ) 1 ' 2 -(d + T)(n 2 - TV 2 )" 2 , (A3) 

where N r =rlf 'is the relative numerical aperture of the rays reaching E in the BFP. Suppose the 
objective illuminated with a Gaussian intensity profile has an effective numerical aperture 
eNA (<NA, see below) corresponding to a radius of the effective pupil r 0 in the BFP plane 
and an effective entrance angle a 0 , the relative numerical aperture could be then normalized as 

N r = r/r 0 eNA = sin a, /sin a 0 eNA . If T=0, then W(r) = {d + Az)(n 2 - N 2 ) 1/2 - d(n] - N 2 ) 1 ' 2 , 

which is the same as [17]. If T=0 and n s =n,, W(r) = Az(n 2 -N 2 r f 2 , as in [36-38]. 

Analytical Zernike coefficients can be computed by expanding W{r) into a series of 
Zernike modes by numerical convolution with the individual Zernike polynomials [21]. The 
AFP is defined as the position where the Zernike defocus is 0. For an objective that meets sine 
condition, all rotationally symmetric aberrations such as defocus, third order spherical 
aberration vary with CG position even when no actual aberrations are present [10]. When the 
CG position is moved with respect to the AFP, Eq. (A. 3) is used to compute the Zernike 
defocus and the Zernike 3rd order spherical aberration as a function of CG position. The 
theoretical slopes of these curves are found to be constant (within 0.2% variation in the ranges 
of depths and CG position studied). 

The aberrations created by moving the CG depend on the effective NA (eNA) of the 
objective illuminated with a Gaussian intensity profile. eNA was defined as the equivalent 
NA of an objective of same focal length illuminated with a constant intensity profile and 
providing the same in-plane resolution. For an objective of pupil diameter D t , illuminated 
with a Gaussian beam of diameter D b (measured at lie 2 ), the truncation ratio t is defined as t = 
D,JD t . Using [31], the 20x/0.5 objective (truncation ration 1.22) has an eNA of 0.47 and the 
63x/0.9 objective (truncation ratio 1.44) has an eNA of 0.85. 
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Appendix B: Supplementary Information 

B.l. CGWS setup 

The components used are: a SLED light source (center wavelength 750 nm, 1.2 mW, FWHM 
bandwidth 23 nm) coupled to a single mode fiber (NA = 0.11, EXS7505-841, Exalos, Swiss) 
with a collimator of focal length 10 mm, AC080-010-B; Thorlabs, US); a polarizer Pol 
(NT47-603; Edmund Optics, US) and a % waveplate (AHWP05M-980; Thorlabs, US); beam 
expander (5x Zoom, BE02-05-B; Thorlabs, US); LI (f = 150 mm, G322352525, Linos, 
Germany), L2 (f = 200 mm, G322353525, Linos, Germany); non-polarizing beam splitters 
BS1 and BS2 (50:50, BS014; Thorlabs, US); reference and sample arm pair of objectives 
(20x/0.5W, UMPlanFI; Olympus) or (63x/0.9 W, HCX APO, Leica); a tilted coverslip 
identical to the one covering the sample to balance dispersion; a Piezo actuator PZT (resonant 
frequency 138 kHz, range 9.1 ± 1.5 (xm, AE0505D08; Thorlabs, US) with a silicon mirror M 
(with low reflectivity in reference arm to improve the contrast of interference pattern) 
mounted on it; motorized linear translation stages TS1 and TS2 (range 28 mm, T-LS28-M; 
Zaber, Canada); a field adjustable diaphragm; a CCD camera CCD2 conjugated to the 
objective focus (resolution 1024x1024, pixel size 14 |xm x 14 um, pixel depth 12 bits, DS11- 
01M15, Dalsa, Canada); a CCD camera CCD1 (resolution 1024* 1024, pixel size 12 |im x 12 
|xm, pixel depth 12 bits, DS-21-01M60-12E, Dalsa, Canada). The pupil apertures of the both 
objectives are conjugated to the active surface of CCD1 by LI and L2. 

The beam dimensions are the following. The collimated beam of the SLED is circularly 
polarized after Pol and QWP. Before entering the interferometer, the beam is expanded to 1 1 
mm (1/e 2 width) by the beam expander, so that the objective pupils are overfilled by a factor 
of 1.22 for the 20^/0.5 objective and 1.92 for the 63^/0.9 objective. The diameter of the 
objective pupils imaged on CCD1 (12 mm for the 20^/0.5 objective and 7.62 mm for the 
63x/0.9 objective) are smaller than the sensitive surface of CCD1 (12.28 mm x 12.28 mm). 

B.2. CGWS Calibration 

We calibrated the CGWS with a real SHS (HASO 3, Imagine Optic) using a mirror as sample 
and with the 20 x /0.5 objective. By displacing the mirror along the microscope axis, we 
introduced a shift of twice the mirror displacement and the defocus aberration at each position 
was recorded both with SHS and CGWS. Both the SHS and CGWS correctly measured the 
slope of Zernike defocus and were in good agreement with the theoretical one (Fig. 9). The 
measured Zernike defocus coefficient of SHS has been normalized in the standard (Noll) 
notation [21]. 

CGWS may report erroneous astigmatism (depending on the size of the scatters) if linearly 
polarized light is used [10], and these errors can be avoided by the use of circularly polarized 
light. Our setup, therefore, contains a quarter wave-plate in the sample arm (Fig. 1) to ensure 
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Fig. 9. CGWS calibration. Measured defocus slope with SHS and CGWS compared to the 
predicted curve (for eNA = 0.47 using the 20x/0.5 objective). 
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that the sample was illuminated by circularly polarized light. We checked in this case that 
spurious astigmatism disappeared, thus demonstrating that no alignment-induced astigmatism 
was present in our system (data not shown). 

B.3. Animal preparation 

All surgical procedures were in accordance with the European Community guidelines on the 
care and use of animals (directive 86/609/CEE, CE official journal L358, 18th December 
1986), French legislation (decret 97/748, 19th October 1987, J. O. Republique francaise, 20th 
October 1987), and the recommendations of the CNRS. 

Before surgery, the rat was anesthetized by urethane injection (1.5 g/kg). Supplementary 
dose of urethane was applied if necessary. The brain was taken out, put into the solution 
(NaCL: 150 mMol/L, KC1: 2,5 mMol/L, Hepes: 10 mMol/L pH = 7.4) and could be stored for 
3 days at most in 4°C refrigerator. For CGWS measurement, the brain was sliced, then held 
on a glass slip, and covered with a cover slip to protect it. During the experiments, the slice 
was always immerged in the solution and the experiment lasted less than 3 hours. 

B.4. Speckle dimension as a function of CG position at shallow depth 

At the smallest depths, the size of the speckle is maximal when the defocus is minimal, as 
expected from the fact that the focus is the smallest at this position (Fig. 10). 




Fig. 10. Speckle size as a function of the CG position relative to the AFP for different depths 
for the 20x/0.5 (a) and the 63x/0.9 (b) objectives. 
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